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Abstract: Ryu and Takayanagi conjectured a formula for the entanglement (von Neumann) 
entropy of an arbitrary spatial region in an arbitrary holographic field theory. The von 
Neumann entropy is a special case of a more general class of entropies called Renyi entropies. 
Using Euclidean gravity, Fursaev computed the entanglement Renyi entropies (EREs) of an 
arbitrary spatial region in an arbitrary holographic field theory, and thereby derived the 
RT formula. We point out, however, that his EREs are incorrect, since his putative saddle 
points do not in fact solve the Einstein equation. We remedy this situation in the case of 
two-dimensional CFTs, considering regions consisting of one or two intervals. For a single 
interval, the EREs are known for a general CFT; we reproduce them using gravity. For 
two intervals, the RT formula predicts a phase transition in the entanglement entropy as a 
function of their separation, and that the mutual information between the intervals vanishes 
for separations larger than the phase transition point. By computing EREs using gravity and 
CFT techniques, we find evidence supporting both predictions. We also find evidence that 
large- symmetric-product theories have the same EREs as holographic ones. 



Contents 



1. Introduction 1 

2. Entanglement Renyi entropy: review 4 

2.1 Basic definitions and properties 4 

2.2 Perturbative MRI 6 

2.3 Replica trick 7 

2.4 Single interval in a CFT 8 

3. Holographic entanglement entropies 9 

3.1 Ryu-Takayanagi formula 9 

3.2 Fursaev's ERE calculation 12 

4. Mutual Renyi information between two intervals 14 

4.1 General properties 14 

4.2 Prediction from Ryu-Takayanagi formula 16 

4.3 Universality in the large-c limit? 18 

4.4 MRI forn = 2 19 

4.5 MRI forn > 2 23 

4.6 Expansion in x 25 

5. Generalizations, open questions, and discussion 28 

A. Symmetric-product CFT computations 30 

A.l Analysis of four-point function of twist operators in C^/Z2 30 

A. 2 Computation of certain OPE coefficients in / Sn 34 



1. Introduction 

The concept of holography originated as an idea about quantum information, that the number 
of qubits that can be stored in a region of space is fundamentally limited by its surface area 
in Planck units. Modern holographic theories go beyond a mere counting of states, and posit 
that the physics governing certain spacetimes can be fully described by a quantum field theory 
residing on its boundary. However, the way that those qubits are organized remains unclear 
on both sides of the correspondence. On one side, we don't yet understand how the states 
are organized in quantum gravity; on the other, despite an in-principle understanding of the 
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state space of quantum field theories, in practice we have to deal with a strongly coupled 
theory with a large number of degrees of freedom. And, of course, the map between the two 
descriptions remains deeply mysterious. 

A useful probe of physical information in quantum systems is the entanglement entropy 
(EE). Here we imagine decomposing a system into two subsystems, A, A'^, with a correspond- 
ing decomposition of the Hilbert space Ti = Ha ® T^A" ■ Given a density matrix p for the 
full system, the reduced density matrix pA, which acts on TIa, is defined by tracing p over 
7iA'= and represents the effective density matrix for an observer who has access only to the 
subsystem A. The EE for A is then the von Neumann entropy of pA- Sa = — tr{pA^i^ Pa)- 
A non-zero EE may be due to the full system being in a mixed state, to information about 
the state being lost by the inability to observe the rest of the system, or to a combination 
of the two effects. The degree of correlation (both classical and quantum) between disjoint 
subsystems may be quantified by their mutual information Ia,b = Sa + Sb — Saub, which 
puts an upper bound on correlators between operators in A and in B [1]. 

In a quantum field theory, it is natural to consider subsystems that are spatial regions. 
Their EEs and mutual informations then tell us about the spatial distribution and correlations 
of quantum information in a given state. Unfortunately, EEs in quantum field theories are 
notoriously difficult to calculate, mainly because one does not have a good way to represent 
the operator InpA- On the other hand, if the density matrix for the full system can be 
represented by a path integral (as in the vacuum or a thermal ensemble, for example), then 
both the reduced density matrix pA and its positive integer powers can also be represented 
in a fairly simple way by path integrals. If those path integrals can be computed explicitly 
for all n, then one can obtain the EE indirectly as follows. Defining the entanglement Renyi 
entropy (ERE) ^J"^ = (lntrp^)/(l - n) for n > 1, one analytically continues in n and 
takes the limit n — )• 1 to obtain the EE. This procedure is called the replica trick. Aside from 
being easier to calculate than the EE, the EREs are of interest in their own right, as a more 

(n) 

refined characterization of the reduced density matrix pA- In fact, knowing S)^ for all n is 
equivalent to knowing the full eigenvalue distribution of pA- In Section 2, we review the basic 
properties of entanglement and Renyi entropies. 

Even given the replica trick, exact results for the EE in field theories are known only in 
very simple cases, such as a single interval in the vacuum of an arbitrary two-dimensional 
conformal field theory [2]. For two disjoint intervals, the EE, and hence mutual information, 
remain unknown even for a theory as simple as that of a compact free scalar [3]. It might 
therefore seem hopeless to dream of knowing the EE in a strongly coupled, large- field 
theory. Remarkably, however, Ryu and Takayanagi (RT) proposed a simple, elegant, and 
universal formula for the EE of an arbitrary spatial region in an arbitrary holographic field 
theory [4, 5]. Their formula, which applies to any state described by a static classical geometry, 
says simply that the EE equals one quarter the area in Planck units of the minimal surface in 
the bulk ending on the boundary of the region A. If correct, the RT formula is not only very 
useful as a calculational tool, but also a significant hint regarding quantum information in 
holographic theories, and probably in quantum gravity more generally (see for example [6]). 
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The RT formula passes several non-trivial checks. For example, it correctly reproduces 
the EE for a single interval in a two-dimensional CFT. A general derivation, using the replica 
trick, was offered by Fursaev [7]. He found that the ERE equaled one-quarter the 
minimal-surface area, independent of n. The analytic continuation in n was thus trivial, 
giving agreement between the resulting value of Sa and the RT formula. In computing 
the ERE, Fursaev performed the necessary path integrals using Euclidean quantum gravity. 
Unfortunately, as we show, the bulk geometries that he used to evaluate the partition function 
are not actually saddle points of the gravitational action. As a result, the ERE he derived 
is incorrect, as we can see by comparing it to the known exact result in the case of a single 
interval in a two-dimensional CFT. We show how the latter result can be reproduced using the 
correct saddle-point action. The RT formula and Fursaev's proof are reviewed and discussed 
in Section 3. 

The question thus arises of whether, in cases where the correct value is not already 
known, we can compute the ERE in a holographic theory, both for its own sake and in order 
to confirm or refute the RT conjecture. Unfortunately, to do so in complete generality, as 
Fursaev attempted, appears to be quite difficult. Therefore, in Section 4, we focus on a 
simple but non-trivial case: two disjoint intervals in a two-dimensional CFT. The RT formula 
predicts a rather interesting phase transition for the mutual information between the two 
intervals as a function of their separation. In particular, for separations larger than a certain 
critical value, the mutual information vanishes, implying a decoupling between the degrees 
of freedom in the two regions. (This behavior of the mutual information is a completely 
general prediction of the RT formula, applying essentially to any two regions in any state of 
any holographic theory. It is closely analogous to the factorization property for disconnected 
Wilson loops [8].) 

The ERE for two disjoint intervals can be expressed in terms of the partition function 
on a certain Riemann surface of genus n — 1. For n = 2, we thus need the torus partition 
function, which fortunately is known for a general holographic CFT [9] . Indeed, as we show, 
it exhibits a phase transition at precisely the same separation as that predicted for the EE by 
the RT formula. For higher values of n, while the partition function is not known explicitly, 
we show using symmetry arguments that the ERE continues to have a phase transition at the 
same separation. This strongly suggests that the same will hold for n = 1, confirming this 
prediction of the RT formula. 

The fact that we can compute the ERE explicitly only for n = 2 precludes analytically 
continuing it to ?i = 1, to directly confirm or refute the full EE predicted by the RT formula. 
We therefore pursue a different strategy. Using the OPE, we expand the ERE, for any given n, 
in powers of the inverse separation between the intervals. The coefficient of any given power 
can be computed explicitly for all n using formulas for conformal blocks, and analytically 
continued to n = 1. We carry this out for a number of coefficients, finding that, thanks to 
a rather intricate pattern of cancellations, in each case the continuation to n = 1 vanishes, 
precisely as predicted by the RT formula. 

As a byproduct of our analysis of the ERE for two disjoint intervals, we find that the result 
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for certain non-holographic CFTs with large central charges, such as large-A^ symmetric- 
product theories, is precisely the same as for holographic ones. It seems that there is some 
form of large-c universality operating here, with a large class of such CFTs having identical 
EREs. This possible feature of the ERE deserves further study. 

We conclude in Section 5 with a list of open questions and possible generalizations of our 
work, and some remarks concerning our current understanding of the RT formula. 

An appendix contains some calculations in symmetric-product theories whose results are 
used in the main text. 

2. Entanglement Renyi entropy: review 

In subsection 2.1 we briefly motivate, define, and state (without proof) the important prop- 
erties of the entanglement Renyi entropy and mutual Renyi information. In subsection 2.2, 
we illustrate these ideas in the simple example of two subsystems that are weakly coupled 
to each other. In subsection 2.3 we then briefly review the replica trick for computing the 
entanglement Renyi entropy, and in subsection 2.4 apply it to the simplest field theory ex- 
ample, a single interval in a two-dimensional conformal field theory. For more details, we 
refer the reader to the books [10, 11] and the review [12]; the latter provides a comprehensive 
introduction to Renyi and entanglement entropies in two-dimensional CFTs. 

2.1 Basic definitions and properties 

Given a density matrix p and a positive real number a 7^ 1, the Renyi entropy is defined as^ 

= ^^Intrp". (2.1) 
1 — a 

At a = 1 the Renyi entropy is defined by taking the limit, and equals the von Neumann 
entropy: 

5 = 5^^) = lim5(") = -tr(pln/)). (2.2) 

Two other interesting limits are limQ_5.o S^°'^ = In dim ^occupied) where ^occupied is the image of 
p, and limQ,_!.oo 5*^"^ = — Inpmax, called the min-entropy, where pmax is the largest eigenvalue 
of p. The following properties of S^'^^ are straightforward to prove: (1) 5'(") > 0, with equality 
if and only if p represents a pure state; (2) S^"'^ is constant if and only if p is proportional to 
the identity on "Hoccupiod) and is otherwise a decreasing function of a; (3) for a > 1 it satisfies 

If the system contains a subsystem A — for example, in a field theory, A could be a 
spatial region^ — then the Hilbert space H can be expressed as the tensor product of Hilbert 

^Writing e-^'°' = {p"-'-)l^'-"~^\ we see that the definition of the Renyi entropy is similar to that of the L'' 
norm of a positive function, ||/||p = (/ fP^^P. The difference is that, whereas in the norm we evaluate the 
integral with respect to a fixed measure, in the Renyi entropy we evaluate it with respect to the very density 
matrix whose entropy we are computing. 

^By region we technically mean co-dimension zero submanifold (possibly with boundary). In this paper we 
will not consider other types of sets, such as single points, fractals, etc. 



- 4 - 



spaces corresponding to A and to its complement A^: H = Ha ® Ha^- Let pA = tr'^^c P 
be the reduced density matrix, defined in TIa, obtained by tracing p over T-Ia'^', this is the 
effective density matrix for an observer who has access only to A. Its Renyi entropy = 
(1 — a)~^ IntiCA Pa is called the entanglement Renyi entropy (ERE) of A, with the special 
case Sa = simply called the entanglement entropy (EE). It can be shown that, if the full 
theory is in a pure state, then = S^^J . 

The EE (but not the ERE for a ^ 1) satisfies an important property called strong 
subadditivity [13, 14], namely, for any two subsystems (or spatial regions) C and D, 

Sc + Sd > Scud + Scnn , Sc + So > Sc\d + -S'dvc • (2-3) 

Strong subadditivity in fact characterizes EE, in the sense that any measure of entangle- 
ment that satisfies (2.3) for all subsystems C and D (as well as certain basic reasonableness 
requirements) must equal the EE [15, 16]. As a special case, strong subadditivity implies 
subadditivity, namely for disjoint subsystems A and B,^ 

\Sa-Sb\<Saub<Sa + Sb. (2.4) 

The second inequality is saturated if and only if the density matrix paub factorizes: paub = 
PA® Pb- Motivated partly by this fact, the mutual information (MI) is defined by 

Ia,b = Sa + Sb - Savjb , (2.5) 

which quantifies the extent to which the degrees of freedom of A and B are correlated with 
each other, including both quantum entanglement and classical correlations. For example, 
the MI puts an upper bound on the correlator between any operator Aa in subsystem A and 
any operator Ab in subsystem B [1]: 

As a consequence of strong subadditivity, the mutual information is monotone under restric- 
tion: if B' C B then Ia,b' < Ia,b- 

The natural generalization of the mutual information defines the mutual Renyi informa- 
tion (MRI): 

Aa) _ q(o) Q{a) nia) „x 
^A,B~'^A "^B '^AUB- y^^-') 

Unlike the MI, the MRI is not necessarily positive. However, it is non-zero only when paub 7^ 
Pa'® Pb, and in this sense still quantifies the extent of correlation between A and B. 

Another reason to study the MRI (including the MI) is that, when we are considering 
a field theory, it is universal, whereas the ERE (including the EE) is cutoff- or regulator- 
dependent. Specifically, when ^ is a spatial region, usually suffers from an ultraviolet 



^The ERE satisfies jS^' — S'^''\ < S^^^ for classical distributions, but not in general for quantum density 



matrices. 
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divergence proportional to the area of the boundary of A. However, if the two regions A and 
B are disjoint and mutually disconnected, then those divergences cancel in the MRI. Since 
the UV regulator generally violates conformal invariance, it follows that in conformal field 
theories the MRI is generally conformally invariant while the ERE is not. Also, in the CFT 
case the ERE suffers from an infrared divergence when one of the regions is infinite in size, 
but this cancels in the MRI (although not when both A and B are infinite). We will see 
explicit examples of these statements throughout this paper. 

In the field-theory context, the restriction monotonicity of the MI implies that the bound 
(2.6) will be strongest when the operators Aa, Ab "cover" the regions A, B respectively. In 
particular, the bound is trivially satisfied for local operators, since then the denominator on 
the left-hand side is generally divergent.^ 

2.2 Perturbative MRI 

Before tackling the computation of entanglement entropies in field theories, as a warm-up we 
first consider the perturbative computation for subsystems that are weakly coupled to each 
other. We will see that in this case the a ^ 1 MRI between the subsystems is parametrically 
larger than the MI, a result that foreshadows the results of Section 4 concerning holographic 
systems. 

We begin by considering a single system, and the effect on its Renyi entropy of a small 
perturbation to its density matrix: 

P = P(0) + V(i) , (2-8) 
where tv P(^o) — P{i) = Oi ^^'^ A is a small parameter. To first order in A we have: 

5'°' = lntrp^„, + ,- ''^yff^ oiX') (2.9) 

1 — a ^ ' 1 — a tr 

S = - tr(p(o) Inp(o)) - Atr(p(i) Inp(o)) + 0{X^) . (2.10) 

Now suppose our system is composed of two subsystems, and the unperturbed density 
matrix factorizes: 

P{o) = PA®PB- (2.11) 

(For example, A and B could represent distinct regions of momentum space, and A the 
coupling constant, in a weakly-coupled field theory.) At zeroth order in A the MRI of course 
vanishes. To first order we have: 

- U + g - « g) ) + 0,A') . (2.12) 



*One could worry about the case where the field theory has a UV cutoff — so the denominator on the left- 
hand side of (2.6) is not strictly infinite — but the MI vanishes exactly. In this case the bound is satisfied for 
a different reason, namely that the density matrix factorizes, so all correlators factorize and the numerator 
vanishes. As far as we know, such a situation can only occur if the space where the field theory lives is 
disconnected. We thank M. van Raamsdonk for pointing out this possibility. 
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For a ~ 1, the operator in the inner parentheses is proportional to {a — 1)^. Hence, to first 
order in A, the MI vanishes: 

Ia,b = 0{X^) . (2.13) 



It can be shown that the order term generically does not vanish. Note that the left-hand 
side of the inequality (2.6) is at most of order A^, so it is reasonable that the right-hand side 
would be of the same order. 



2.3 Replica trick 

Unfortunately, in practice there are very few known methods for computing EREs (or EEs) 
in field theories. One of the most useful is the replica trick, which we will review below, 
that allows one to compute the ERE 5^"^ for integer n > 1 [2]. In favorable circumstances 
a simple analytic form for 5^"^ for general real a can be found which fits those data points, 
and from this form the EE can be read off by setting a = 1.^ It is important to say at the 
outset that in proceeding this way we are merely presuming to have guessed the ERE S^°'^ 
correctly; firstly, nothing guarantees (in an infinite-dimensional Hilbert space) that the ERE 
is analytic, and, secondly, the values of a function on a countably infinite set (in this case, 
the integers larger than 1) are not sufficient to fix a unique analytic continuation. (There 
exist analytic functions, such as (1 — a)~^sin7ra, that vanish for all integer a > 1 but not 
elsewhere, including at a = 1.) Having stated this caveat, for the rest of the paper we will 
assume that all EREs we consider are indeed analytic functions of a. We will find nothing 
inconsistent with this assumption. 

The replica trick applies when the theory is in a state, such as the vacuum or a thermal 
state, whose partition function can be obtained by a path integral over some Euclidean 
spacetime E (possibly with some operator insertions, which for the purposes of this discussion 
we will consider to be part of E). Let ^ be a spatial region, and En the n-sheeted cover of 
E with the sheets connected along branch cuts placed at ^ on a constant Euclidean-time 
surface. Then tiAP^ ~ ^n/Zf, where Zn is the partition function of the theory on En (and. 



^Throughout this paper a wiU he in the interval [0, oo], while n will be a positive integer. 
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in particular, Z\ is the partition function for the original theory).^ Hence we have (for n > 1) 

To be more concrete, let us further specialize to a two-dimensional conformal field theory 
and let A be the union of disjoint intervals where Ui < Vi < Uj+i. Then we 

can rewrite the expression (2.14) in terms of correlators of twist operators in the symmetric- 
product theory C^/Sn (or equivalently C"/Z„), computed on E: 

= j^ln{<jl{ui)aU{vi) . . .al{uN)<jU{vN)) • (2.15) 
1 — n 

The correlator of twist operators is divergent, due to the singular geometry of En at the 
branch point. It can be regularized by regularizing each twist operator separately; hence the 
notation (t| and o"!^, where e is a UV cutoff length, for the regularized twist operators [17]. 

2.4 Single interval in a CFT 

As an example of the application of (2.15), the ERE for a single interval, in the vacuum, is 

where c is the central charge of C and Cn is a scheme-dependent quantity. Here we used the 
fact that the twist operators have scaling dimension 

d,= i(„-i). (2.17) 

The (simplest) analytic continuation of (2.16) to non-integer a is^ 

^{a) C ( . 1\. ( V - U 



^If the theory contains fermions then one needs to specify their boundary conditions across the constant 
Euclidean-time surface where the sheets are sewn together, which we will call S. The original partition function 
Z\ is computed with a sign-flip on the fermionic flelds across 5*. More generally, Zn is computed with a flip on 
SV\A^ on each sheet of iJ„, along with one on S n A where the nth sheet connects to the first sheet (but not 
on the other n — 1 copies of 5* PI A). Hence a curve that winds around a branch point n times, returning to the 
same point on crosses n + \ sign flips. The resulting overall flip for even n is part of the deflnition of the 
twist operators g\ and cr_i of the next paragraph. When is evaluated by passing to a coordinate system 
that is single- valued on En, this overall flip is canceled by the branch cut in the coordinate transformation for 
the fermionic fleld. (For example, in complex coordinates if 2 is a local coordinate on E with the branch point 
at 2 = 0, and t — z^''" is a single-valued local coordinate on E„, then ^pt = {dz/dt)^^^^pz = n^^^'^t^^~"^^^iJz- 
For even n the factor t'""!)/^ has a branch cut with a sign flip.) Hence there is no operator insertion in the 
new coordinate system. However, there may still be sign flips around non-contractible cycles. (See for example 
the case of the torus in footnote 14 below.) 

^AU CFTs will be assumed compact, unitary, and modular-invariant in this paper. 

*It is interesting to ask what eigenvalue distribution for p[u,v] gives rise to the a-dependence S^^\-^ = 
(1 -I- l/a)C seen in (2.20) (neglecting the subleading and scheme-dependent quantity Ca)- This question can 
be answered by deflning a fictional "Hamiltonian" H = — lnp[„ „] acting on Hiu^v]- Then S'^^\ is related to 
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which yields the EE [2] 

5[„,,] = ^hi(")+ci. (2.21) 

Note that S^l^^j indeed satisfies the properties (1), (2), (3) mentioned below equation (2.2). 

It is also possible to obtain the result (2.16) (and thereby derive the scaling dimension 
(2.17)) by computing Z„ and applying (2.14). The computation of Z„, is done as follows [17]. 
We are in the vacuum, so the Euclidean spacetime E is simply the plane, to which we add a 
point at infinity to make it topologically a sphere. The multi-sheeted surface En is then also 
topologically a sphere. A Weyl transformation maps the metric ds^ on En to a fiducial metric 
ds"^ = e~'^ds^ on the sphere. We then have Zn = e^^Z, where Z is the partition function of 
C on the sphere with the fiducial metric, and 5^ is the Liouville action: 

SL = ^j g^'^ {9'"'d^<pd,<P + 2i?0) (2.22) 

(which depends on C only through its central charge) . For n > 1 the metric on En has a conical 
singularity at each branch point u,v, so the Liouville action is divergent. The divergence can 
be regulated by replacing a disc of radius e about each branch point with a smooth metric, 
which defines the regularized twist operators cr^i- 



3. Holographic entanglement entropies 
3.1 Ryu-Takayanagi formula 

In this subsection we will provide a brief summary of Ryu and Takayanagi's proposal for the 
entanglement entropy (EE) in field theories with holographic duals [4, 5], along with some of 
the evidence supporting it. A more complete review can be found in [19]. 

The Ryu-Takayanagi (RT) conjecture is a proposed formula for the EE of a given spatial 
region A in certain states of holographic field theories whose dual gravitational theory is 
classical Einstein gravity (possibly with matter). Specifically, the proposal concerns states 

the free energy of H at the temperature a~^: 

F = -i mtre-'^- Sr -{^-^)^- (2.18) 

(The first two equahties apply to the Renyi entropy of any system.) The density of states that gives rise to this 
temperature dependence for the free energy is easily found, in the saddle-point approximation, by performing 
a Legendre transform: 

^(^)^{!!:p(2cv^(i^-c)v^), ^>c' ^'-''^ 

where E is the eigenvalue of 6. (See [18] or [12] for the form of the full inverse Laplace transform.) It is 
interesting that (up to a shift of E by C) p{E) has the same form as the Cardy formula for the asymptotic 
density of states in a CFT on a circle. Note that we have not determined which physical observable H 
represents — it is not necessarily related to a physical energy. 
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that admit a description as static classical solutions in the dual theory, such as the vacuum 
and thermal states.^ We work in a fixed constant-time (i.e. timelike-Killing-field orthogonal) 
slice of the bulk. The conjecture states that 

area(mA) , . 

where rriA is the minimal-area surface in the bulk that is homologous to A, i.e. such that there 
exists a region va with dr^ = AU rriA- (We use the term "homologous" loosely, given that 
A is not generally closed. This topological condition was originally suggested by Fursaev [7], 
and plays a crucial role in several checks of the proposal.) The area is evaluated with respect 
to the Einstein-frame metric. 

An interesting question, assuming the RT formula is valid, is how it gets corrected by 
quantum effects and by higher-derivative (e.g. a') corrections to the classical action in the 
bulk. Quantum effects presumably lead to Gn corrections to (3.1) (starting at order G^), 
although a specific form has not been proposed. In the presence of higher-derivative correc- 
tions to the classical bulk action, it is expected, based on consistency with black-hole entropy 
(discussed below), that in the formula (3.1) the area is replaced by Wald's black-hole entropy 
formula [21] (or at least some functional that coincides with it on horizons). 

The RT proposal passes several basic checks. For example, if A is the entire boundary, 
then Sa should simply be the statistical entropy of the state. Indeed, according to the RT 
proposal we should take ruA to be the minimal surface in the bulk that is homologous to 
the boundary; this will generally be the horizon, if there is one, giving agreement with the 
Bekenstein-Hawking entropy. If there is no horizon, then the boundary is homologically trivial 
in the bulk (i.e. the topological boundary of the bulk is precisely the boundary where the 
field theory lives); hence the minimal surface is the empty set, giving Sa = 0.^*^ (Again, this 
is the order entropy — the RT formula does not capture the entropy due for example to 
a gas of gravitons, which is of order G^.) Furthermore, when the total entropy is zero (or of 
order G^), then if we instead take A to be a subset of the boundary, we expect from (2.4) 
that Sa = Sa<^- Indeed, in this case the entire boundary is homologically trivial in the bulk, 
so A and A'^ are homologous, implying = m^c 

Another important check on the RT proposal is that it satisfies the strong subadditivity 
(SSA) property (2.3) for any regions G and D, as can be shown by a simple geometrical 
argument [23]. (Interestingly, the proof of SSA based on the RT formula is far simpler 



^Possible generalizations to time-dependent states were proposed in [20]. 

^"Since we are considering static spacetimes, any black holes in the spacetime should be eternal, so the max- 
imally extended spacetime may include other, topologically disconnected boundaries. Consider, for example, 
the maximally extended spacetime of the AdS-Schwarzschild black hole, whose boundary has two connected 
components. It is believed [22] that the field theory defined on both boundaries represents the thermofield 
double of the field theory defined on only one boundary. In this picture, the black hole spacetime, which 
represents a thermal and therefore mixed state in the original field theory, represents a pure state in the ther- 
mofield double. This result is faithfully reproduced by the RT prescription; the full boundary (including both 
components) is homologically trivial in the bulk, giving Sa = 0. 
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than the general proof.) Since, as mentioned in subsection 2.1, SSA characterizes the EE, 
this is quite strong evidence in favor of the RT formula. (However, it is not sufficient to 
prove its correctness, since it only shows that (2.3) is satisfied for subsystems corresponding 
to geometrical regions, whereas for the characterization proof one needs it to hold for all 
subsystems.) The proof extends trivially to the inclusion of higher-derivative corrections, as 
long as they are extensive. 

As a final check, let us see how the RT formula reproduces the EE (2.21) of a single interval 
[u,v], in the vacuum of a two-dimensional CFT. The vacuum is described holographically by 
AdSa, whose metric on a constant-time slice is 

ds^ = §^Uz' + dy^). (3.2) 

here y is the coordinate along the boundary and z is the radial coordinate, with the boundary 
being at z = 0. We employ a simple UV cutoff in which we shift the boundary to z = e. 
The minimal surface m[u,v] is a geodesic connecting the points on the boundary (y, z) = 
(n, e), (i), e), which is an arc of a circle (almost a semi-circle) with center ((n + i))/2,0). 
Applying (3.1) and using the standard holographic relation t-Ads/G^ = 2c/3, one finds [4, 5] 

matching (2.21). (In this scheme, the finite part ci vanishes.) In higher dimensional CFTs, 
although one does not have exact formulas for the EEs even of simple regions, the leading 
UV divergence is known and matches that predicted by the RT formula [4, 5]. 

For a time it was believed that the RT formula should only apply to connected regions. 
(See for example the paper [24].) The reason was that, when applied to the union of two 
intervals (a case that will be considered in detail in the next section), it disagreed with a 
calculation by Calabrese and Cardy [25] which (like the formula (2.21) for a single interval) 
was supposed to be valid in any two-dimensional CFT. However, those same authors have 
since shown that their original calculation was incorrect. At present, there is no reason to 
believe that the RT formula, if it is valid at all, would not apply equally well to connected 
and to disconnected regions. For example, all the checks discussed above apply to both cases 
(including the last one, which can be considered a computation of the EE of the disconnected 
region (— oo,u] U [f,oo)). 

When applied to a disconnected region, the RT formula makes a fascinating prediction 
for the mutual information (MI) between its components, similar to the phase transition for 
disconnected Wilson loops found by Gross and Ooguri [8]. For simplicity, let us consider 
two disjoint and mutually disconnected regions A, B. Each has a corresponding minimal 
surface m^i, niB and region r^i, r^. (We assume the generic situation that ta^tb are disjoint 
and mutually disconnected.) When we consider the region Au B, the disconnected surface 
rnA^ruB is topologically allowed and locally minimal. Assuming that the full bulk spacetime 
is itself connected, surfaces will also exist that connect A and B. However, if the separation 
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between A and B is sufficiently large compared to their sizes (and any other scales defined 
in the theory or state), then rriA Dtub will necessarily be the globally minimal surface. Then 
we have Saub = Sa + Sb, so the MI Ia,b vanishes. More precisely, Ia,b is of order G^, 
rather than Gj^^.^^ This implies that, from a quantum information point of view, the two 
regions are approximately decoupled from each other. (See for example the bound (2.6) on 
correlators between A and B. Note however that this bound does not give us new information 
about correlators of local operators, for which the left-hand side is always of order G^.) This 
prediction is striking in its generality, applying as it does to all holographic theories. (In 
Section 4 we will find evidence that in fact it applies even more generally, to a large class of 
large- field theories.) If we then imagine bringing A and B closer to each other, then it 
may happen that, at some critical separation, the minimal surface will switch from uia Ums 
to one that connects dA and dB (see for example figure 1). In this case, the MI will (in 
the thermodynamic/classical limit Gn — 0) undergo a first-order phase transition; it will 
become non-zero, with a continuous value but discontinuous first derivative as a function 
of the separation between A and B. Section 4 will be devoted to a detailed study of these 
phenomena in the simplest example, namely two intervals in the vacuum of a two-dimensional 
CFT. 

3.2 Fursaev's ERE calculation 

In the paper [7], Fursaev gave a derivation, based on the replica trick, of the RT formula. 
In this subsection, we will briefly summarize his argument, and then point out a flaw that 
results in an incorrect value for the entanglement Renyi entropy (ERE). 

In our sketch of Fursaev's argument, for simplicity we will take the bulk action to be pure 
Einstein gravity; matter fields and higher-derivative (e.g. a') corrections are straightforwardly 
incorporated, as he discusses. We will also assume that the ultraviolet divergence in the 
field theory is cut off in some manner whose details will not concern us. Fursaev's starting 
point is (2.14), where Zn is the partition function on the n-sheeted Euclidean spacetime 
Recall that, if A is the spatial region whose EE we are computing, then the sheets of En 
are connected by a branch cut along A on a constant-time slice. In a holographic theory, 
this partition function is given by the gravitational path integral over Euclidean geometries 
whose conformal boundary is En- In the classical limit, this path integral goes over to its 
saddle-point approximation e~'^™'", where Smin is the minimal value of the Euclidean Einstein- 
Hilbert action among extrema obeying the boundary conditions. Fursaev constructs a set of 
geometries with boundary En, then minimizes the Euclidean action within that set. He takes 
as given the bulk Euclidean spacetime E representing the original state of the system; its 
boundary is E and its Euclidean action is — InZi. He takes n copies of E and connects 
them along a branch cut r^, which is a spatial region in E lying in the same constant-time 
slice as A. In order for this n-sheeted bulk geometry to have boundary En, the part of the 

^'^A closely related phenomenon, in which the A and B are held fixed but the bulk spacetime is deformed, 
was discussed in the paper [6]. 
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boundary of rA that lies in E must coincide with A (i.e. drA H -E = A); apart from this 
condition, the choice of is at this point arbitrary. The branch "point" is niA, the rest of 
the boundary of {rriA = Ova \ A and drA = Au ruA)- He now evaluates the Euclidean 
Einstein-Hilbert action for this geometry. There are two contributions. First, the geometry 
is made up of n copies of E, so there is a contribution —nlnZi, which is independent of 
the choice of ta- In addition, the Ricci scalar has a delta function along the branch "point" 
niA, which is co-dimension 2 and hosts a conical singularity with excess angle 27r(n — 1). It 
therefore contributes a term (n — 1) area(myi)/(4GN) to the action. Minimizing this action 
over all possible choices of rA, he obtains the minimal surface niA, and (from (2.14)) the ERE 

(„) _ area(mA) . . 

- 4Gn ■ 

Since there is no n-dependence, the analytic continuation is particularly simple: = 
area(myi)/4GN- Finally, setting a = 1, he obtains the RT formula (3.1). 

The problem with this derivation is that the action has been extremized only with respect 
to a subset of the degrees of freedom in the metric, namely the choice of r^. The resulting 
field configuration is therefore not guaranteed to be a true saddle-point, and in fact it does 
not solve the Einstein equation: the Einstein tensor has a delta function supported on ruA 
due to the conical singularity, with no corresponding source. 

We can confirm that the ERE (3.4) is incorrect by comparing it to the exact result in a 
case where the latter is known. For example, when j4 is a single interval in the vacuum of a 
two-dimensional CFT, the exact ERE (2.16) depends on n (by the factor 1 -|- 1/n), whereas 
the Fursaev result (3.4) is independent of n. What is the true saddle point in this case? The 
Euclidean space £" is a plane, and the corresponding bulk geometry E is hyperbolic 3-space 

(a.k.a. Euclidean AdSs). The saddle point corresponding to the n-sheeted cover En is also 

The easiest way to see this is to add a point at infinity to E to make it a sphere; then its 
n-sheeted cover £"„ is also a sphere, so the corresponding bulk geometry is H"^. This geometry 
is smooth, in contrast to Fursaev's, which is n copies of glued together in such a way as to 
create a conical singularity along the geodesic connecting the endpoints of A. Given that the 
bulk geometry is for all n, why does its action depend on n? The bulk action is divergent 
due to the infinite volume near the boundary; while the full bulk geometry is for any n, 
the cutoff geometry depends on n. (The full geometry depends only on the Weyl class of the 
boundary metric, which is the same for all n, since the sphere admits a unique Weyl class. 
On the other hand, the cutoff geometry is sensitive to the actual boundary metric. This is 
the holographic manifestation of the Weyl anomaly [26].) The n-dependence can most easily 
be calculated by performing a Weyl transformation to put the metric on En into a standard 
form and taking into the account the resulting change in the partition function due to the 
Liouville action, as described at the end of subsection 2.4, or equivalently by the holographic 
renormalization procedure [26]. 

It is worth noting that the true saddle point with boundary En can be obtained 
topologically by gluing n copies oi E = together in precisely the manner described by 
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Fursaev. This will continue to be the case in the more complicated examples we will study 
in the next section, suggesting that, while it carries the wrong metric, Fursaev's construction 
may be topologically correct in general. This would explain why the topological condition on 
the minimal surface uia that he suggested — that rriA should be homologous to A — appears 
to be correct. 

Finally, it is intriguing that, while Fursaev's value (3.4) for the ERE is incorrect for 
n > 1, it somehow manages to give the right answer for the EE (n = 1), assuming that 
the RT conjecture holds. We can only speculate that, if there is some sense in which the 
spacetimes En and their bulk duals can be defined for non-integer values of n, then his 
construction may be correct "at linear order" in a neighborhood of n = 1. 

4. Mutual Renyi information between two intervals 

As we saw in subsection 2.4, the entanglement entropy for a single interval in the vacuum 
of a two-dimensional CFT depends only on the theory's central charge. The fact that the 
Ryu-Takayanagi formula correctly reproduces this entropy, as shown in subsection 3.1, is 
an important check on the proposal, but does not give us any new information. The next 
simplest configuration we can consider in such a theory consists of two disjoint intervals. As 
suggested by the fact that the Renyi entropies (2.15) depend in this case on four-point rather 
than two-point functions of twist operators, we would expect the EE to depend on the full 
operator content of the theory, rather than simply its central charge. As we will see in this 
section, the RT formula can give us significant new physical information in this case. The 
new predictions in turn give us the opportunity to subject the formula to new and highly 
non-trivial quantitative tests. 

We begin by reviewing the necessary formulas and setting up the basic properties of the 
ERE for two intervals. 

4.1 General properties 

We consider two separated intervals [ui,t;i], [^2,^2] ("i < vi < U2 < V2) in the vacuum of a 
conformal field theory C with central charge c. As discussed in subsection 2.1, it is convenient 
to consider the mutual Renyi information (MRI) between the two intervals, 

[ui,Dl],[u2,t'2] ["lif^l] [«2,«^2] [«l,fl]U[n2,ll2] ' \ ■ ) 

which measures the extent to which the degrees of freedom of the two intervals are entangled 
with each other (including both classical correlations and quantum entanglement). 

We first consider the integer case a = n > 1. Using (2.15), the MRI is given in terms of 
a finite ratio of four-point and two-point functions in the symmetric-product theory / Sn'- 

(„) -^inf /";(";^"-)(";;"^("^^"-^(y^^M (4.2) 

[u^,v{\,[u2M n-1 \{al{ui)a'^-^{vi)){al{u2)(Jt,^iv2)) J ^ ' 

1 / {(Ti{ui)a^i{vi)ai{u2)a-i{v2)) \ ^-^ 



n-1 \{ai{ui)a_i{vi)){ai{u2)(T-i{v2)) 
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where we've defined the renormahzed twist operators: 

" W(0).t'(i)>v. • 

This is an example of the UV divergences in the EREs, which occur at the endpoints of the 
intervals, cancelling in the MRI, as discussed at the end of subsection 2.1. 

Since the twist operators are primaries, the transformation law for the four- and two- 
point function implies that I[ui,vi],[u2,v2] is conformally invariant, and therefore depends only 
on the cross-ratio 

. ^ , (4.5) 

{U2 - Ui){V2 - Vi) 

which lies in the interval < x < 1. By a conformal transformation, the four points ui, vi, 
U2, V2 can be brought to 0, x, 1, oo respectively, so we have: 



^^"^(^) - Cd.oo] = (^''" (ai(0)cT_i(x)ai(l)a^,(oo)>) , (4.6) 



where a'_i{oo) = lim^_>oo z^'^" a-i{z). (The scaling dimensions d^j of the twist operators are 
given by (2.17).) Notice that, like the UV divergence, the IR divergence in the ERE cancels 
in the MRI. The four-point function, and therefore I^"'\x), is an analytic function of x in the 
interval < x < 1. 

It is useful to note that the four-point function in (4.6) can be expanded as a power series 
in X, where the powers are the dimensions dm of operators A'^ in the orbifold theory,^^ and 
the coefficients are given in terms of OPE coefficients: 

\ m 



n 



Note that only untwisted operators contribute to the sum. Assuming we are dealing with a 
unitary theory, the operator with lowest scaling dimension is the unit operator, for which (by 
the normalization of the twist operators) the OPE coefficients are 1. Hence /(")(x) goes to 
as X — )• 0, as we would expect on physical grounds. For example, if we fix the sizes of the 
intervals and take their separation to infinity, we would expect all correlations between them 
to go to zero. We will study the higher-order terms in the expansion (4.7) in subsections 4.4 
and 4.6, and in the appendix. 

A final important property of the MRI, implied by the invariance of the four-point func- 
tion in (4.6) under x — >• 1 — x, is 

/(")(1 -x) = l(")(x) + + In ^—^ . (4.8) 



6 V I X 



^^Throughout the paper, we use primes on (untwisted) operators of C" /Sn, to distinguish them from oper- 
ators of C. 
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At the level of the definition (4.1) of the MRI, this relation is due to the fact that, in a pure 



state (in this case, the vacuum), S 



(a) 



4"\so s^r^ 



s, 



(a) 



(a) 



'A '^A'' ' '^^ ^[0,x]U[l,oo] '^[-oo,0]U[a:,l] ^[0,l-x]U[l,oo] ' 

We have listed five general properties that the MRI satisfies for integer a > 1, but for 
the reasons given we either know or expect each to hold for general values of a: 

1. UV finiteness, and IR finiteness when one of the intervals is semi-infinite; 

2. conformal invariance, implying 



i: 



{a) 

[ui,Vl],[u2,V2] 



[O,x],[l,oo] 



(4.9) 



(where < x < 1); 



3. 



lim l(")(x) = 0; 



4. for all X, 



l(")(l_a;) =/(°)(x) + ^ fl + -^ 

o \ a J 



In 



X 



(4.10) 
(4.11) 



5. analyticity of I^°'\x) function of x. 



So far we have not assumed anything about the theory C (other than unitary and com- 
pactness). In the rest of this section, we will study the function I^°'\x) in holographic CFTs, 
as well as certain other theories with large central charge. 

4.2 Prediction from Ryu-Takayanagi formula 

As in the holographic derivation of the EE for a single interval, reviewed in subsection 3.1, 
we use the fact that the holographic dual of the vacuum is AdSa, with ^Ads/^N = 2c/3, and 
we cut off integrals near the boundary at radial coordinate value z = e. 




Figure 1: The two locally minimal surfaces for the boundary region U [u2,'y2]- The global 

minimum is the one on the left is when x < 1/2, and the one on right when x > 1/2, where x is the 
cross-ratio defined in (4.5). 

The RT formula is straightforward to apply to the union of two intervals [ui, vi] U [u2, ^2]- 
There are two locally minimal surfaces in the bulk that are homologous to this boundary 



- 16 - 



region, as shown in figure 1. The first is the union of the minimal surfaces for the two 
intervals separately, m^is = m^ui.vi] U ''t^[u2,v2] (similarly for the corresponding bulk region 
^'dis = U r^u2,v2])- This has "area" (i.e. length) 

area(mdis) = area(m[„j^^j]) + area(m[„2^^2]) (4.12) 

= 2^Aas In ( ^1 + 2i^,s In ( ^1 (4.13) 



(see (3.3)). The other locally minimal surface connects ui to V2 and U2 to vi: rricon = 
"^["1,^2] ''^[vi,u2]- (The corresponding bulk region is a semi-annulus connecting the two 
intervals: rcon = i~iui,v2] \ ''^lvi,u2]-) Its area is 

area(mcon) = 21^43 In (^"^^^^ + 21^63 In (^''l^^ . (4.14) 

It is easy to see that mjis is the globally minimal surface when x < 1/2, and mcon otherwise 
(see (4.5)), so 

Sluuvi]ulu2,v2] = ^^niin(area(mdis),area(mcon)) 

^ iln{{vi-ui){v2-U2)/e^), x < 1/2 
3 \ln{{v2-ui){u2-vi)/e^), x > 1/2 ' 

Combining (4.15) with (3.3), we obtain the following mutual information: 

(..^ ^ fo, X < 1/2 ^ ^ 

[ (c/3) ln(x/(l — xj) , X > 1/2 

Of the five properties of the MI listed at the end of the previous subsection, this formula obeys 
the first four. It does not obey the last — analyticity — as it has a discontinuous first derivative 
at X = 1/2 and vanishes for x < 1/2. These two features were anticipated in the discussion 
at the end of subsection 3.1. The discontinuity in the first derivative occurs because the 
global minimum switches between the local minima as we vary x, and is reminiscent of phase 
transitions due to competing saddle points of the Euclidean action, such as the Hawking-Page 
transition. As in that case, the transition is presumably sharp only in the classical limit in 
the bulk, which corresponds to the thermodynamic (c — )• oo) limit of the CFT, and gets 
smoothed out by finite-c effects. Similarly, the vanishing of the MI for x < 1/2 is presumably 
true only at order c; if the MI vanished exactly for x < 1/2, then the reduced density 
matrix for the two intervals would factorize, P[ui,vi]u[u2,v2] — P{ux,vx\ ® P\u2,v2\-: implying that 
the two intervals are completely decoupled from each other, a rather unphysical situation 
(in particular, it would violate the inequality (2.6), as we will discuss below). Thus we 
should expect both perturbative and non-perturbative corrections to (4.16) in Gn ~ c~^, 
with the first perturbative correction at order (F' . Nonetheless, since the MI is apparently 
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parametrically small for x < 1/2 — smaller than the EE for either interval separately or for 
their union, and smaller than the MI for x > 1/2 — it appears that the density matrix factorizes 
approximately. 

Unlike quantum corrections, we do not expect higher-derivative (e.g. a') corrections to 
the classical bulk action to change the result (4.16), for the following reason. As discussed 
in subsection 3.1, such corrections are believed to correct the area functional appearing in 
the RT formula without changing the basic prescription of minimizing over topologically 
allowed surfaces. The symmetries of AdSa guarantee that the minimal surfaces shown in 
figure (1) remain uncorrected; furthermore, the corrected "area" of each curve is unchanged 
when written as a function of c, since we know that the EE is always given by (2.21). In fact, 
this argument applies for any bulk gravitational action, not just Einstein- Hilbert with small 
higher-derivative corrections. 

It is interesting to consider the result (4.16) in the context of the general inequality 
(2.6). Local operators in holographic field theories are dual to bulk fields. Smearing such 
operators over the two intervals will give a result of order — no matter how the operators are 
normalized, the numerator and denominator will be of the same order in c. This is consistent 
with the result (4.16) (and shows that it must be corrected at order c^). However, for x > 1/2, 
the right-hand side of (2.6) is of order c. It would be quite interesting to find examples of 
operators ^[^^ ,;^], A[u2,v2] that saturate (2.6), or at least for which the left-hand side is of order 
c. It is then interesting to ask what happens to these operators for x < 1/2. Do they simply 
cease to exist? Or does their two-point function become parametrically smaller? Finding 
explicit examples of operators with such behavior would not just constitute strong evidence 
in favor of the RT formula, but also clarify in what sense the density matrix factorizes. 

4.3 Universality in the large-c limit? 

Given any family of CFTs C that admit a large-c limit, such as holographic ones, we can 
consider the expansion of the MRI in powers of c~^. Since the number of degrees of freedom 
is of order c, the leading term will be at most of that order, so we have 



In particular, we focus our attention on the leading function (x). In the previous subsec- 
tion, we used the RT formula to compute, for holographic theories. 



and argued that this result should hold no matter what the bulk gravitational theory is. As 
discussed, (4.18) has two striking qualitative features, namely its discontinuous first derivative 
at X = 1/2 and the fact that it vanishes for x < 1/2. In the rest of this section, we will 
study if ^(x) using the replica trick, and find independent evidence for both phenomena. 

(2) 

In the next subsection, we will compute (x) in holographic theories, and show that the 



l(")(x) = /i"^(x)c + 4"^(x) + 0(c-^) . 



(4.17) 




X < 1/2 
X > 1/2 



(4.18) 
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result is independent of the details of the bulk theory (e.g. the presence of higher-curvature 
corrections), and applies also to symmetric- product theories C = Cq /Sn, even though their 
large-c limit is not described by classical gravity. Like (4.18), the result will have a phase 
transition at x = 1/2. In subsection 4.5 we will argue that this phase transition occurs also in 
l["\x) for n > 2, at least in holographic theories. Then in the last subsection we will use CFT 
techniques to study the expansion of l[°'\x) in powers of x for general a, and find evidence 
that every coefficient in this expansion goes to in the limit a — t- 1. That analysis will 
assume very little about the CFT C, essentially just that the number of operators below any 
given dimension stays finite as c goes to infinity, a condition that holds for both holographic 
and symmetric-product theories (but not, for example, in the power theory Cq without the 
orbifold) . 

These results not only give strong quantitative support to the RT formula, but point 
to a broader picture, namely that a large class of large-c CFTs — including holographic and 
symmetric-product theories — share the same leading MRI l[''\x), function of both a 
and X. Although we do not know the explicit form of this function except for a = 1, 2, we can 
deduce that it is analytic in x except at x = 1/2, where it has a discontinuous first derivative, 
and satisfies the following properties: 

lim/i"^(x) = 0, (4.19) 

x—^O 

li")(l -x) = /i")(x) + l(l + -)ln — . (4.20) 

o \ a J X 

Based on these considerations, it appears that in the range < x < 1/2 the MRI is 
parametrically larger for a ^ 1 (where it is of order c) than for a = 1 (where it is of order 1). 
This is similar to what we found in the perturbative calculation of subsection 2.2. It would 
be interesting to find a simple toy model of a system with degrees of freedom, in which 
the MRI between two subsystems is of order N, but the MI is only of order 1. 



4.4 MRI for n = 2 

In this subsection we will begin by expressing the mutual Renyi information I^'^\x) (some- 
times called the mutual collision information) in a general CFT in terms of its torus partition 
function. Using this expression, we will calculate the order-c part /| (x) in a general holo- 
graphic CFT, finding that — like the RT prediction (4.18) for l[^\x) — it is analytic except 

(2) 

at X = 1/2, where it has a discontinuous first derivative. We will then show that /| (x) is 
precisely the same function in large-A symmetric-product theories, supporting the idea of 
universality (i.e. theory-independence in the large-c limit) proposed in the previous subsec- 
tion. 

We begin by applying the formula (4.6) for n = 2. In the C'^ jZ^ orbifold theory, there 
is a unique twist operator a = a\= a-\. Lunin and Mathur [17] showed that its four-point 
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function is given by 
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{a{0)a{x)a{l)a{oo)) = (2®x(l - x)) 



-c/12 



Zil , 



(4.21) 



where Zu is the partition function for C on a flat rectangular torus^^ with modular parameter 
T = il; X and / are related by 



X 



eliii) ■ 

As X goes from to 1, I goes from oo to 0, with x 
2). Since 



(4.22) 

1/2 corresponding to Z = 1 (see figure 



1 



(4.23) 



the invariance of the four-point function (4.21), and hence of the ERE, can be traced to the 
modular invariance of the torus partition function, Zu = Z^/i. The reason for the appearance 
of the torus partition function of C is that the four-point function of twist operators is the 
(renormalized) zero-point function on the two-sheeted Riemann surface E2 with a branch 
cut connecting to x and another one connecting 1 to 00, which is a torus with complex 
structure t = il. The Weyl transformation that flattens it, when plugged into the Liouville 
action, leads to the prefactor (2*x(l — x))'^/^"^ . 

Plugging (4.21) into (4.6), we obtain / 



/(2)(x) = lnZi;--^ln 



12 



28(1 -x) 



x^ 



In Zil 



In 



ol(ii) 
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The first term in (4.24) is (— / times) the free energy 
of C on a circle of unit circumference at temperature 
l""^. (Note however that the basic cycles of this torus, 
which we interpret as space and Euclidean time di- 
rections when we speak of the free energy, are not 
the same as the space and time directions of the Eu- 
clidean plane E where the theory was originally de- 
fined, whose double cover is £'2- Rather, the spatial 
circle of the torus encircles the points and x, staying 
on one sheet, while its Euclidean time circle encircles 
X and 1, crossing each branch cut once.) 




Figure 2: Modular parameter t — il 
for the two-sheeted Riemann surface with 
branch points at 0, a;, 1, 00. The rela- 
tion between x and I is given by equation 
(4.22). 



^^In terms of Lunin and Mathur's variables, x — l/io and il = t — — l/rLunin-Mathur- The four-point function 
of twist fields was computed in [27] in the case where the underlying CFT is a free scalar field. 

^''The fermion sign fiips explained in footnote 6 imply that fermionic fields should have antiperiodic (NS) 
boundary conditions on both cycles of the torus. The reason is that, on the double-sheeted plane, in going 
around either cycle one encounters two sign fiips, so no overall flip. The coordinate transformation to the 
torus introduces a flip, just as when passing from the plane to the cylinder. Hence the partition function is 
not invariant under the full modular group, but it is invariant under r — ^ — I/t- 
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-2-Kdl 



+ 




(x < 1) 



(4.29) 



(4.30) 



Let us consider the expansion of (4.24) for small values of x, where x w 16e~'^'.^^ The 
behavior of the torus partition function is universal in this limit (for compact unitary CFTs), 
InZji ~ 27rcZ/12, which precisely cancels the leading behavior of the second term in (4.24), 
giving a vanishing MRI as expected (equation (4.10)). The leading x-dependence depends on 
the gap in the operator spectrum of C. If the lowest non-unit operator A has dimension d 
and multiplicity m (where fermionic operators are counted negatively), then 

27rc/ 
IT 

so 

(2c + m)e~2'^' ~ 2-S(2c + m)x^ 
^^-2iTdl ^ m{x / 16)'^'^ , 

This term can be matched onto the leading term in the expansion in intermediate states (4.7), 
by noting that the lowest-dimension operator of C^/Z2 appearing in the aa OPE, other than 
the unit operator, is A<^ A if d < 1, and the stress tensor if d > 1. The OPE coefficients 
are computed and matched to (4.30) in Appendix A.l (see also the discussion around (4.39)), 
where we also consider more generally the matching between the Lunin-Mathur formula for 
the four-point function (4.21) and its expansion in intermediate states. 

A simple example of the application of (4.24) 
is to a free scalar compactified on a circle of ra- 
dius R. The torus partition function is 



so [28, 3] 
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This is plotted against x for several values of R in 
figure 3. The small- a; behavior is as predicted by 
(4.30), with the lowest non-unit operator having 
dimension d = 1/{2R^) (for R^ > 1) and mul- 
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We now turn to holographic CFTs, briefly reviewing Maldacena and Strominger's result 
for the torus partition function [9]. Expanding the free energy in powers of ~ Gn, the 
leading term is of order c and is given by the Euclidean action of the dominant saddle point. 
Here the boundary condition is simply that the conformal boundary should be the torus with 
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the expansions of the theta functions for smah q are as follows: 
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T = il; there are no operators inserted in the path integral so no fields other than the metric 
are sourced. For / > 1 (x < 1/2) the dominant saddle point is the Euclidean BTZ black hole, 
which is topologically a solid torus in which the Euclidean time circle (the circle of length 
/) is contractible. For / < 1 (x > 1/2) the dominant saddle point is Euclidean AdSa with 
the Euclidean time direction periodically identified; the topology is a solid torus in which the 
spatial circle is contractible.^^ The phase transition between the two saddles, the Hawking- 

(2) 

Page transition, is first-order, so the free energy, and hence (x), has a discontinuous first 
derivative. Specifically, the Euclidean actions of the two saddle points yield [9] 



l2^cVl2 + 0(cO), l>l 



so 



which is plotted in figure 4. Note that, although I\ (x) does not vanish in the region x < 
1/2, it is numerically quite small — smaller than I^'^\x) for the free scalar by two orders of 
magnitude or more. The expansion for small x is 2~^x^; comparing to (4.30), it is as if the 
lowest-dimension operator of C^/Z2 is the stress tensor. In fact, there are other operators, 

(2) 

but since their multiplicity is finite in the limit c — t- oo, they do not contribute to . We 
will discuss this expansion in detail in subsection 4.6. 

If we compare the bulk saddle-point geometries used to derive (4.33) to the ones obtained 
from Fursaev's construction, we see that they are topologically identical but metrically dif- 
ferent. To describe the geometry obtained from Fursaev's construction, we add a Euclidean 
time direction, coming out of the page, to each diagram in figure 1, and consider the double 
cover of the resulting three-dimensional geometry, branched over rdis and rcon respectively. 
On the left-hand diagram, relevant when x < 1/2, the cycle on the boundary that encircles 
ui and fi, staying on one sheet, is contractible through the bulk; this is the spatial circle 
of the torus. On the right-hand diagram, relevant when x > 1/2, the boundary cycle that 
encircles vi and U2, crossing both branch cuts, is contractible through the bulk; this is the 
time circle of the torus. Hence in each case the topology is precisely the same as that of 
the true saddle-point geometry. However, their metrics are different; in particular, while the 
former are singular, the latter are smooth. 

The formula (4.33) for the torus partition function, and therefore the formula (4.34) 
for the MRI, applies not only to holographic CFTs but also to symmetric-product theories 
C = Cq /Sn at large A^, where Co is any (compact unitary) GET [29]. The basic reason is that, 
for I < 1, the effective temperature for the long strings, which dominate the partition function, 
is enhanced by a factor of A^, so the theory is effectively always in the high-temperature limit; 



^^Note that the contractibihty of the two cycles of the boundary torus requires anti-periodicity of fermions 
on both. As explained in footnote 14, these are precisely the boundary conditions we have in this case. 
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Figure 4: Coefficient of c in tlie mutual Renyi information of intervals [0,x] and [l,oo] in a general 
holographic CFT. The two plots differ only in the scale of the vertical axis. In particular, the plot on 

(2) 

the right shows that ![ ' is non-zero (although quite small) for a; < 1/2. 

the partition function for Z > 1 is then given by modular invariance. An example is the 
supersymmetric (T^)^ /Sn theory, which is conjectured to be connected in a moduli space 
to type IIB string theory on AdSa x x T^. It would seem reasonable then to guess that 
the torus partition function is given by (4.33) for all theories on this moduli space. In other 
words, the MRI appears to enjoy a non-renormalization theorem. 

The phase transition at x = 1/2 can be understood in terms of the expansion Zn = 
e~27r(di+c/i2)/^ where the di are the scaling dimensions of the operators of C. In any fixed 
theory, with finite c, this expansion converges and is analytic for all x G [0, 1). In the large-c 
limit, the operators of C can (roughly speaking) be divided into those with scaling dimensions 
of order 1 ( "short strings" ) and those with scaling dimensions of order c ( "long strings" ) . Each 
long-string operator makes a contribution to the sum that is exponentially suppressed in c. 
However, the number of long-string operators is exponentially large in c, so they may actually 
dominate the sum. In fact, whether the short strings or the long strings dominate depends on 
the value of /, and hence of x. In both holographic and large- symmetric-product theories, 
short strings dominate for x < 1/2 and long strings for x > 1/2. Thus the order-c part of the 
free energy is due entirely to short strings for x < 1/2 and to long strings for x > 1/2. 

4.5 MRI for n > 2 

In this subsection we will extend our study of /^"^(x) to larger values of n. Although we 
will not be able to give explicit formulas, we will argue that all the main qualitative features 
carry through. In particular, the existence for all n > 1 of a discontinuous first derivative in 
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/}" (x) at X = 1/2 constitutes significant evidence in favor of the RT formula, which predicts 
precisely such a phase transition for n = 1. 

According to equation (4.6), the MRI /^"^(x) for general n is given in terms of the four- 
point function of twist operators in the symmetric-product theory C"'/Sn- This four-point 
function is in turn equal to the (renormalized) zero-point function of C on the surface 
made of n sheets connected by a branch cut running from to x and another from 1 to cxd. 
This is the Riemann surface for the equation = z(z — l)/(z — x), which has genus n — 1 
and a complex structure that depends on x. As in the genus- 1 case studied in the previous 
subsection, the surface En can be taken by a Weyl transformation to a fiduciary metric with 
the same complex structure, for example the constant-curvature metric. Hence there is an 
analogue of (4.6), in which /^"^ is written as a sum two terms: 

/W(x) = ln4")+c/;j_,,Jx). (4.35) 

(n) 

Zx is the partition function of C on the surface carrying the fiduciary metric. The second 
term is derived from the Liouville action for the Weyl transformation from En to the fiduciary 
metric; aside from the overall coefficient c, it is independent of the particular theory C, giving 

(n) 

a universal contribution to . 

(n) 

The geometrical term I| gcoj^cti-ic has not been explicitly computed for n > 2. However, 
assuming that the fiduciary metric is chosen to depend smoothly on x (as does, for example, 
the constant-curvature metric), the Weyl transformation and hence -^I'g'comctric ^^^^ smooth 
functions of x. 

Meanwhile, the genus-(n — 1) partition function Zx"'^ is known explicitly for n > 2 only in 
a small number of CFTs. For holographic theories, despite considerable progress (especially 
in the context of pure gravity theories), explicit formulas are not available, even in the large- 
c limit; see for example [30, 31, 32, 33, 34]. Even in the absence of an explicit formula, 
however, we can argue that the partition function is smooth except at x = 1/2, where it has 
a discontinuous first derivative, just as we saw for n = 2 in the last subsection. It is known 
that phase transitions, analogous to the Hawking-Page transition, occur at fixed points of 
the mapping-class group (the group of large diffeomorphisms of the Riemann surface). The 
reason is that different saddle points of the bulk gravitational action are effectively mapped 
onto each other by the action of the mapping-class group, and therefore at a fixed point they 
necessarily have the same action. On the surface there is an element of the mapping- 
class group that effectively takes x to 1 — x; for example, the surfaces En with x = 1/3 and 
X = 2/3 have the same complex structure up to the action of an element of the mapping-class 
group. That element permutes the cycle that encircles and x with the one that encircles 
X and 1. It has a unique fixed point, namely x = 1/2. Hence we expect a phase transition 
in (the order-c part of) lnZx"'\ and therefore in l[^\x), at x = 1/2, and only there. This 
is precisely the property predicted by the RT formula for l[^\x) (see (4.18)). (Furthermore, 
it seems likely that, as we saw for n = 2, the dominant saddle has the same topology as 

^^We thank A. Maloney for helpful discussions on this point. 
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predicted by Fursaev, i.e for x < 1/2 the cycle encircling and x is contractible in the bulk, 
while for x > 1/2 the cycle encircling x and 1 is contractible in the bulk.) 

In the previous subsection we used the fact that the torus partition function is the same 
(at leading order in c) for large- symmetric-product theories {C^ /S^) as for holographic 

(2) 

ones, to support the claim of universality (theory-independence) for /) (x). It would be 

(n) 

interesting to investigate whether the same holds for Zx for n > 2. Our results in the next 
subsection, which apply for all n, support such universality. 

4.6 Expansion in x 

(2) 

In subsection 4.4 we derived an explicit expression for /| (x) that applied to both holographic 
and large- symmetric-product theories. Unfortunately, as we discussed in the last subsec- 
tion, the computation of higher-genus partition functions, and therefore l[^\x), remains out 
of reach technically in such theories. Even if we could find explicit expressions, their analytic 
continuation to general a, and in particular to a = 1, may not be feasible. (For example, 
Calabrese, Cardy, and Tonni were able to compute /("'^(x) for n > 1 for a compactified free 
boson [3]. However, the analytic continuation of the resulting expression, which involves a 
Riemann-Siegel theta function, is unknown.) In this subsection, therefore, we will take a 
different approach, and compute l["'\x) order by order in x for general n > 1. The coefficient 
of each power of x will be a simple enough function of n to allow us to analytically continue 
it straightforwardly.^^ As discussed in the previous subsection, for each n > 1, l[^\x) is 
analytic on [0,1/2], but not on larger intervals. Assuming that this property continues to 
hold for l[°'\x) for general a, the formulas we derive will be valid on that interval. We are 
particularly interested in testing two hypotheses concerning the coefficient of each term in 
the power-series expansion of /j"^(x): that it is "universal" in the sense of subsection 4.3, i.e. 
the same for all theories in the class we are considering; and that it goes to zero in the limit 
a —7- 1, in agreement with the RT prediction (4.18). We will find significant evidence in favor 
of both hypotheses. 

We are considering theories C for which the number of operators with any given dimension 
is finite in the large-c limit. This condition applies to holographic and Cq / Sn theories (where 
Co is held fixed as we take c ~ A^ — )■ oo), but not for example to the theory Cq without the 
orbifold. We will also assume that the ?i-point functions of primaries in C do not diverge in 
the large-c limit; again, this property holds for holographic and Cq /Sn theories. 

We begin with the relation (4.7), in the form 

exp ((n - l)/(")(x)) =Y,c^\,mc"'a,a^,x''- . (4.36) 

m 

At this stage we have not taken the large-c limit, and the convergence of the OPE dictates that 

the sum on the right-hand side converges and is analytic for all x G [0, 1). Only untwisted- 

sector operators occur in the aia-i OPE; these are of the form A'^ = {Ai-^ (8) • • • (8> ^j„)sym, 

^*We remind the reader that throughout this paper we use n to denote a positive integer and a a non-negative 
real number. We assume that all quantities are analytic functions of a. 
^^This procedure was applied to the compactified free boson in [3]. 
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where the Ai are operators of C ("sym" means average over aU permutations; if some of the 
Ai are identical, then this is the same as the average over distinct permutations). 

Before considering the large-c hmit, it is interesting to ask what the assumed analyticity 
in a of I^'^\x) impHes about the coefficients on the right-hand side of (4.36). First, except 
for the leading term, 1, the total coefficient of each power of x must be a multiple of n — 1 
(i.e. when analytically continued must have a root at n = 1). Although we don't know a CFT 
proof of this statement, it can be tested to any given order. For example, for the holomorphic 
and antiholomorphic parts of the stress tensor, T' and T' , we have 

C^V.T' = c-^^^f, = y . (4.37) 

The Zamolodchikov metric on these operators is Qt'T' = Gfif/ = nc/2, since the central 
charge of the orbifold theory is nc; they do not mix with each other or with other operators, 
so g^'T' =gf'f' ^2/{nc), and 

C <j\(J—\ — C cri(J-_i — ■ l4.oo) 

nc 

Hence their contribution to the sum on the right-hand side of (4.36) is 2aa:^, where 

dl in' - l?c 



2nc 288n3 



(4.39) 



which has the required factor oin—1. The contribution of the operator {Oi®Oi®I®- ■ ■(^l)syrm 
where Oi is a primary of C, is computed in Appendix A. 2 (for integer values of the scaling 
dimension of Oi); the result again has a zero at n = 1. More generally, analyticity of 
I^"\x) demands that the contribution of an operator made up of k non-unit operators of C, 
{Ai^ <8) • • • (8) Ai,. (8) / (8> • • • ^ I)sym, should contain a factor (n — 1) (n — 2) • • • (n — A; + 1) , simply 
because that operator only exists for n > k. 

We now consider the large-c limit. Since we are working order by order in x, we will 
consider only "short-string" operators of C, i.e. those whose scaling dimensions are finite in 
the large-c limit. A convenient machinery for systematically computing terms in the sum 
(4.36) is provided by the conformal blocks. We thus write it as a sum over primaries O'^ (of 

exp ((n- l)/(")(x)) =Y,Cm.Hhm,x)J^Chm,x)x''^+''"^ , (4.40) 

■m 

where Cm = c'^VimC™crio-_i , hm and hm are the weights of O^, and T{hm, x) is (up to a factor 
of x^'^""^'^) the conformal block with all 4 external operators of weight h^j = hcj = da/2. {Cm 
and the conformal blocks also depend implicitly on c and n.) The conformal blocks can be 
computed straightforwardly (albeit somewhat tediously) to any desired order in x. Using 
Mathematica'^^ , we computed them to order x^. In the following expression, the first two 



For these computations we used the package Virasoro.nb, available at 
http : / /people .brandeis. edu / " headrick / physics / . 
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terms are exact, while the subsequent ones have been expanded in powers of a ^ ~ c 
F{hm,x) = 1 



+ [a + f2Q + 0{a-^)) 



+ ( (l + ^)a + /3o + 0(a-^))x3 
+ ( + /4i« + ho + 0{a-^) ) 



+ 0{x'^). (4.41) 



The fij are rational functions of and n (regular at n = 1); their precise form is not 
important for us, except for one feature we will point out below. We now factor out the 
positive powers of a (i.e. of c), since those determine l["'\x). It turns out that they organize 
themselves naturally into an exponential: 

T{hm, x) = F{hm, x) exp (ax^ + ax^ + g4ax^ + g^ax^ + 0{x^)) , (4.42) 

where, by definition, F contains only non-positive powers of a: 

F{hm,x) = l + ^x+ (/20 + O(a-i)) x2 + (/30 + 0{a'')) x^ + 0(x4) . (4.43) 

Remarkably, thanks to some cancellations among the fij, the coefficients of ax^ and ax^ turn 
out to be independent of hm- 

km 1309n4 - 2n2 - 11 
54 - /4 1 - /20 - ^ 

hm 589n^ - 2n^ - 11 
95 = hi - ho - J20 - 54— = T20n^ ' ^ ' 

Assuming that this pattern continues to higher orders, it allows us to pull the exponential 
out of the sum (4.40), and write: 

/f)(x) = J(«)(x) + ^""'llfj^^' {x' + x'+ 94x' + g^x' + 0{x')) , (4.45) 

where j("^(x) is the contribution to l[^\x) (if any) from the OPE coefficients Cm- 

(x) = lim ^ In ( V CmF{hm, x)F{hm, x)x'^'"+'^- ) . (4.46) 

n — 1 c^-oo c \^ I 

Before discussing j(")(x), let us point out several noteworthy features of the second term 
of (4.45). First, it does not depend at all on the particular theory, supporting the universality 
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proposed in subsection 4.3; this is a consequence of the cancellation of the /im-dependence in 
54 and §5, (4.44). Second, if we set n = 2, it agrees with the expansion to fifth order of (4.34); 
hence J^^^ (x) vanishes at least to fifth order. Third, it can be straightforwardly continued to 
non-integer values of a, and vanishes at a = 1: 

l[^\x) = jW(x) + ^" {x' + + g^x^ + g5X^ + 0{x^)) . (4.47) 

The fact that the second term vanishes at q = 1, which provides strong quantitative evidence 
in favor of the RT formula, can be traced to the fact that the conformal block (4.41) depends 
on c through a ~ (n — l)^c. 

It remains to ask what we can say about J('^)(x). Since we are disallowing theories in 
which the number of primaries of a given dimension is proportional to c, J^") {x) will be non- 
zero if and only if some of the coefficients Cm contain positive powers of c. Some (but not 
all) of the primaries of / Sn are products of primaries of C: O'^ = (Oii ■ ■ ■ C'j„)sym- For 
these, as we show in Appendix A. 2, the OPE coefficient c^'^aim is given by the (symmetrized) 
n-point function in C of Oj^ , . . . , Oi„ , where the operators are placed at distinct nth roots of 
unity (see (A. 19)). In holographic and / S]^ theories, these n-point functions go like c^~^f'^ , 
where k is the number of non-identity operators, so Cm ~ (?~^ (this is for k > \; for /c = 0, 
i.e. the identity of /Sn, Ci = 1, while for k = 1, Cm = 0, since the one-point function of a 
non-identity operator vanishes). For example, for k = 2 we have a two-point function, which 
is clearly independent of c {Cm. is computed in this case in appendix A. 2). Hence primaries 
of C^/Sn that are products of primaries of C do not contribute to J^"'\x). However, there are 
other primaries of C^/Sn that are made up of descendants of C. The simplest example is 

(^l + ^^L_ia, . (4.48) 

The computation of Cm for such operators is more involved, because of the more complicated 
transformation law of the C-descendants {L-iOi and L-iL-iOi, in the above example) in 
going from the n-sheeted plane to the standard plane. It is possible that such operators 
contribute to J^'^\x). If so, the conjecture is then that J^"\x) is independent of the particular 
theory C, and vanishes at q = 1. It should be straightforward to compute Cm and test these 
conjectures in specific examples, such as (4.48). We leave this to future work. 

In this subsection, we have provided non-trivial evidence, based on the expansion (4.36), 
that l["'\x) is theory-independent and, for < x < 1/2, vanishes at a = 1. In view of the 
pattern we have found, it would clearly be desirable to have some general understanding of 
the structure of the OPE coefficients at large c that leads to these properties. We leave the 
exploration of this structure to future work. 

5. Generalizations, open questions, and discussion 

In the previous section, through the study of Renyi entropies, we provided strong evidence 
in favor of the Ryu-Takayanagi formula. We focused on one of the simplest non-trivial field- 
theory examples, namely two disjoint intervals in the vacuum of a two-dimensional GET. 
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Along the way, we found evidence that a large class of large-c theories share the same en- 
tanglement (Renyi) entropies. It would be interesting to extend our analysis to more general 
situations, including: more than two intervals; states other than the vacuum, such as thermal 
states; CFTs on the circle rather than the line; CFTs in more than two dimensions; and 
non-conformal field theories. In particular, it is clear that the two key qualitative predictions 
of the RT formula persist in all these examples, namely that there is a phase transition in the 
mutual information between two regions as a function of their sizes and separations, and that 
it vanishes on one side of the phase transition. One should be able to test these predictions 
using similar techniques to the ones used in this paper, namely classical gravity and the OPE. 
One should also be able to test whether the EREs are the same for non-holographic theories 
with large central charges. 

Our analysis leaves a number of open questions. We will start with the more technical 
ones, and move towards the more conceptual. First, our calculation of l["'\x) in subsection 
4.6 left out the term J^"'\x), which comes from primary operators of C^/Sn, such as (4.48), 
that are composed of descendants of C. It would be useful to evaluate this term, at least up 
to some power of x, to confirm both its theory-independence and that it vanishes at n = 1. 
More generally, it should be possible to understand on general GET grounds the pattern 
found in subsection 4.6 that, in the four-point function of twist operators, every factor of c is 
accompanied by a factor of (n — 1)^. 

Second, it would be very interesting to compute the MRI l["'\x) explicitly for n > 2 in 
holographic and large- symmetric-product theories, to see, first, if they agree, and second, if 
they indeed have a phase transition at x = 1/2. Better yet would be to analytically continue 
the resulting expressions to general a, and directly confirm or refute the RT formula in this 
case. 

Third, as discussed in subsection 4.2, the RT formula formula predicts that the mutual 
information between the two intervals is of order c when they are close together, but only 
of order 1 when they are far apart. It would be quite illuminating to find operators A[ui,vi]i 
A[u2,v2] that saturate the inequality for x > 1/2 (2.6) (or at least such that the left-hand side 
is of order c), and to understand what happens to them for x < 1/2. A related question is 
raised by the fact that, for x < 1/2, the MRI is of order c (for a ^ 1). In order to understand 
the behavior better, it would be interesting to find a simple toy-model system with a large 
number of degrees of freedom, in which the MRI between two subsystems is of order of the 
number of degrees of freedom while the MI is only of order 1 . 

Lastly, and perhaps most importantly, we should ask what the status of the RT formula 
is, given the results of this paper. On the one hand, we have provided strong evidence that 
it is correct. On the other hand, have we understood any better why it should be true? In 
particular, why does the minimal surface play a critical role in the entanglement entropy, and 
what is the physical significance of the bulk region that it bounds? Eursaev's proof [7], 
though incorrect, had the advantage of explaining in a simple and elegant manner the role of 
the minimal surface. On the other hand, in the Renyi entropy calculations we have performed 
in this paper, this role is not so clear. Rather, the agreement between the Renyi entropies 
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and the RT formula appeared to be almost fortuitous. Clearly, while the RT formula provides 
a tantalizing hint about the structure of quantum information in holographic theories, most 
of that structure still remains hidden from view. 
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A. Symmetric-product CFT computations 

A.l Analysis of four-point function of twist operators in C^/Z2 

In this appendix, we consider a general modular-invariant, compact, unitary CFT C with 
central charge c, and its symmetric product C^/Z2. (At the end we also make some comments 
about the general symmetric product C"/S'„.) The orbifold theory has central charge 2c, and 
contains a single twist operator a with conformal weights h„ = h^j = c/16. Lunin and Mathur 
[17] computed the four-point function of these operators, showing that it is determined by 
the partition function Zr of C on a flat torus with modular parameter r, 

{a{0)a{x)a{l)a'{oo)) = \2^x{l - x)\~"^'^^ , (A.l) 

where r and x are related by^^ 

The reason for the appearance of the torus partition function of C is that the four-point func- 
tion of twist operators is the (renormalized) zero-point function on the two-sheeted Riemann 
surface E2 with a branch cut connecting to x and another one connecting 1 to 00. E2 
is a torus with complex structure r. The Weyl transformation that flattens it leads to the 
prefactor |2'^x(l - 

Both sides of (A.l) can be decomposed into a sum of states, and we would like to 
understand the relationship between these two decompositions. The torus partition function 

^^In terms of Lunin and Mathur's variables, x — 1/w and r = — l/rLunin-Mathur- In the bulk of the paper 
we consider x to be real and lying in the interval < x- < 1, but in this appendix we will let a; be a general 
complex number. As in the main text, a'i{cs3) = limz^oo z"^" ai{z). 
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is a sum over states Am in C: 



Zr = ^g'^™-/24^-^™-e/24 = ^ ^^ ,^(q)^^ j^^(q) , (A.3) 
m i 

where q = e^'^*'^. In the second equahty, we have grouped the states into conformal famihes. 
Each family is labelled by its primary operator Oi, and Xc,hi is its Virasoro character: 

oo 
N=0 

where d{N) is the number of descendants of Oi at level A^. The decomposition (A.3) can be 
obtained by cutting the torus along a cycle and inserting a complete set of states. In the 
usual presentation of the torus as C/(Z + rZ), that cycle should be horizontal. 

Meanwhile, the left-hand side of (A.l) can be written as a sum over intermediate states 
A'l of C2/Z2, with weights {h[, h[): 

(a(0)a(x)a(l)a'(oo)> = ^ c,,ic' ,„x^'^~''^x^'^-''^ , (A.5) 

I 

where Ccrai = (f'''(oo)f''(l)-^;(0)) ^^"^ '^'o'o" coefficient of A'l in the a-a OPE. Assuming 

for clarity that |x| < 1, this decomposition is obtained by cutting the sphere on a circle of 
radius r {\x\ < r < 1) around the origin, which separates the twist operators located at and 
X from those located at 1 and oo, and inserting a complete set of states. 

The intermediate states in (A.5) can also be organized into conformal families, leading 
to a sum of conformal blocks. However, since the conformal families of C^/Z2 are not in 
one-to-one correspondence with the conformal families of C, and we are trying to reproduce 
the sum (A.3) which is over the latter, we will organize the intermediate states slightly 
differently. First we note that only untwisted states appear in the sum, and these are of the 
form A'l = Am "X" An + An ^ Am, where Am, An are states of C. We are inserting this state 
on the circle of radius r mentioned in the previous paragraph. In the presence of the twist 
operators, we can consider that we are working in the theory C on the Riemann surface E2, 
where the circle is two circles, one on each sheet; we are inserting An on one circle and Am on 
the other. These two circles both represent the same cycle of the torus, namely the horizontal 
cycle mentioned below (A. 4). In other words we have cut the torus into two finite cylinders. 
Each cylinder has An inserted on one boundary and Am inserted on the other. 

We now perform the Weyl transformation that turns E2 into the flat torus. Two things 
will happen. First, we get the geometrical factor |2^x(l — x)\~'^/^'^, as computed by Lunin 
and Mathur, which is independent of the states. Second, each state gets mapped by the 
action of the conformal group to a linear combination of states. By definition, this group acts 
within conformal families. Hence if An and Am are not in the same family, then the cylinder 



For simplicity we are taking all states to be bosonic. 



amplitude vanishes. So we can gather the terms in (A. 5) into conformal families of C: 

{x,x), (A. 6) 

i 

where 

Am , An descendants of Oi 

(The set of operators in C^/Z2 of the form Am ® An + An ® Am where Am and An are both 
descendants of the primary Oi in C, is the miion of several conformal families of C^/Z2. Hence 
Ki includes several conformal blocks of C^/Z2.) Each term of (A. 6) corresponds to precisely 
one term in the sum on the right-hand side of (A. 3), and the Lunin-Mathur formula tells us 
that 

K,{x,x) = |28x(1-x)|"'/''xca('?)Xca(^)- (A-8) 

It is interesting that Ki factorizes as a holomorphic times an antiholomorphic function. 

Each state in the sum (A. 7) contributes to Ki a monomial in x, x, while each state in the 
sum (A. 4) contributes to Xc,hiiQ) a monomial in q. The complicated mixing between states 
due to the action of the conformal group is reflected in the complicated relationship between 
X and q. However, the leading terms for small x on the two sides of (A. 8) can be matched 
easily. On the right-hand side the leading term is due to the primary Oi itself, so we have 

where we used the expansion for small x, g ~ 2~^x^. The leading term on the left-hand 
side is due to the operator O ■ = Oj ig) Oi, which has weights {h'^,h[) = {2hi,2hi). 0[ is 
primary, so (taking it to be normalized in the Zamolodchikov metric) we have Co-o-j' = c^' aa = 
((j(0)C'-(l)(j'(oo)). To evaluate this three-point function, we consider the theory C on the 
two-sheeted Riemann surface with a branch cut extending from to oo, and with Oi inserted 
at the point z = 1 on both sheets. We can use the map z = io relate this to the two-point 
function {Oi{—l)Oi{l)) in the t-frame. (The factor arising from the Weyl transformation is 
absorbed in the renormalization of the twist fields.) All in all we find 

W = = 2-4''-4^^ (A.IO) 

which leads immediately to agreement with (A. 9). 

In principle equation (A. 8) can be checked to higher orders. Consider, for example, the 
conformal family of the identity. For convenience, let us divide both sides of (A. 8) by the 
leading term: 

\x\^/^Ki{x,x) = 



x^ 



-c/24 



Xc,oiq) 



(A.ll) 
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oo ^ 



Generically, the conformal family of the identity is a fuh Verma module except the states 
L_i|0) and L_i|0) and their would-be descendants, which vanish. (For the minimal models 
there are also other missing states.) In that case the character is 

°° 1 1/24 /I \ 

so the holomorphic part of (A. 11) is 

n=2 

The first few terms in the expansion in powers of x are: 

1 + 2-^cx^ + 2~^cx^ + 2~^^(c2 + 465c + 2)x^ . (A.14) 

In the expansion in states of C^/Z2, (A. 7), the quadratic term is due to the stress tensor, 
while the cubic term is due to L_3|0) = dT. The correct matching of the coefficient for the 
former can be seen by setting n = 2 in (4.39). 

Modular invariance means that the torus partition function can be written as a sum of 
characters in a different way, namely 

i 

where q = e~^'^*/^. This decomposition is produced by cutting the torus along its "vertical" 
cycle. Meanwhile, associativity of the OPE means that the four-point function of twist 
operators can be decomposed in intermediate states with each state contributing a power of 
1 — x (instead of x as in (A.5)), by cutting along a circle centered on 1 that separates 1 and 
X from and oo. That circle corresponds to two circles on £'2, both representing the vertical 
cycle. Thus the two decompositions can be mapped to each other just as we did above. It 
is interesting that the associativity of the OPE in C^/Z2 is directly related to the modular 
invariance of C. 

If we attempt to generalize the above analysis to the analogous four-point function of 
twist operators 

(c7i(0)(j_i(x)ai(l)c7'„i(oo)> (A.16) 

in the / Sn orbifold theory, the following structure emerges. The Riemann surface £"„, has 
n sheets joined by a branch cut extending from to x and another one extending from 1 
to 00. This surface has genus n — 1, and the circle centered on 0, that separates the twist 
operators located at and x from those located at 1 and 00 decomposes into n circles, which 
separate £"„ into two n-punctured spheres. Again, only untwisted states, which are of the 
form {Am-^ (^■■■® Am„)sym; enter in the sum we insert on that circle. We are left with a sum 
of squares of n-point functions of C (to be contrasted with (A.5), which is a sum of squares 
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of three-point functions of C^/Z2, or in this case C^/Sn)- Unhke in the C^/Z2 case, the Ann 
do not all have to belong to the same conformal family of C to contribute to this sum. For 
this reason, this decomposition is less immediately useful than for the case n = 2. There will 
also be an overall geometrical factor coming from the appropriate Weyl transformation. 

A. 2 Computation of certain OPE coefficients in C"/S'„ 

In this appendix we will consider primary operators in the symmetric-product theory C^/Sn, 
of the form 

where the are primaries of C. We will first show that c'^'^aim, the OPE coefficient with the 
twist operator ui, is given in terms of the n-point function (in C) of the component operators 
Oi. We will then focus on the simplest non-trivial case, with only two non-identity operators 
(necessarily the same, for c°"Vim to be non-zero): 

0'^ = iO,(^Oi---(^I(^---(^I)^^^ . (A.18) 

We calculate Cm = c'^^aimC^aia-i in the case where the scaling dimension di of Oi is an 
integer, showing that its analytic continuation to non-integer n vanishes in the limit n — )• 1. 
We begin with the more general operator (A.17). We compute: 

aim ^cT—imfJi 

= (aii(0)aj(l)>^:/^^ (cTli(0)O:„(l)af (oo)>^„/^^ 

■^^'''^ ((0,,{e^-'h0^,{e'-'/n■■■0^M)) ) . (A.19) 

^ \ / C/ sym 



'C on K 



sym 



n 



(All correlators except the one marked "C on En" are evaluated on the Riemann sphere.) In 
the first line we used the fact that the twist operators are normalized, and both they and O'^ 
are primary. In the fourth we used the definition of the twist operators to move to the original 
theory C on the n-sheeted surface En, where the sheets are connected by a branch cut running 
from to oo (the positions of the twist operators). The operator Oi^ is positioned at 1 on the 
jth sheet, denoted e^'^*-'. In the last line we conformally mapped En to the plane by t = z^/". 
The geometrical factor from the associated Weyl tranformation is independent of the operator 
insertions, cancelling the factor (cj1^(0)(t|(1))~^. The operator positions are mapped to the 
nth roots of unity. Since each Oi - is primary, under the conformal transformation it becomes, 
in the t-frame, \dz/dt\ '^^^Oi. = n '^'^Oi- (we are assuming for simplicity that the operators 
are scalars; including spin there are additional phase factors, which depend on the position 
and therefore have to be included before the symmetrization over permutations). 
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To compute Cm one also needs to know the Zamolodchikov metric for O'^- A bit of 
combinatorics shows that this is 

ni!---nfc! 

ymm = j , (A.2Uj 

n! 

where k is the number of types of operators Oi in O!^ and the n; are the number of operators 
of each type (so Yli ^^i = n). 

We now specialize to the operator (A. 18). From (A. 19) we obtain 



^ ' j.*.— 1 



^ n—l 

n - 1 ^ I 



-2di 

(A.21) 



We wish to analytically continue this expression in n.^^ We were not able to do this for 
general dimension di, but in the next paragraph we will show that, when di is a (positive) 
integer, the sum in (A.21) is a polynomial in n of degree 2dj, with a root at n = 1. Hence 
c'^Vim is a rational function of n that is regular at n = 1. Since Gmm = 2/{n{n — 1)), 
Cm = ^rnm(c'^^o-im)^ a rational function with a zero at n = 1. This was the statement that 
was used in subsection 4.6. 

In order to analytically continue the sum in (A.21), we note that the summand equals 
the reside of the pole at t = e^'^*-'/" of the function 

f{i) = —, -r-^ ■ (A.22) 

We are assuming that di is a positive integer, so /(t) is single- valued. It is easy to show that 
it is regular everywhere on the Riemann sphere except for a pole at each nth root of unity. 
In particular, at t = 1 there is a pole of order 2dj + 1, and the sum in (A.21) equals minus its 
residue. Writing u = t — 1, this is the coefficient of n^'^*"^ in the expansion of the function 

Now, it is clear that this coefficient is zero for n = 1, since u~^{\ + u)*^'"^ has no term of 
order v?'^^~^ . It remains to show that it is a polynomial of degree 2di. To do this we re- write 
(A. 23) as 

When we expand the sum in large parentheses in powers of u, the leading term is u, and 
after that the coefficient of n*" is a polynomial in n of degree m — 1. It follows that, when 



^^The analytic continuation of the sum in (A.21) was also considered in [3]. In particular, an expression was 
derived that allowed numerical approximations to be computed. 
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we expand the whole expression in powers of u, the leading term is u^^, and after that the 
coefficient of u"^ is a polynomial in n of degree m + 1. So in particular the coefficient of v?''^^~^ 
is a polynomial of degree 2di . 
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